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1 Introduction 

Radiative corrections to jet production in e + e~ annihilation were computed a long 
time ago QTJ, 0, 0. These calculations were, however, performed for massless quarks. 
In most practical applications this is sufficient, since, at relatively low energy, the b 
fraction is strongly suppressed, and at high energy (i.e. on the Z peak and beyond) 
mass effects are suppressed. Nevertheless there are several reasons why a next-to- 
leading-order calculation is desirable. First of all, at sufficiently high energies, top 
pairs will be produced and mass effects there are very likely to be important. A second 
reason is to understand the relevance of mass corrections, due to bottom production, 
to the determination of a s from event shape variables. As a third point, quantities 
such as the heavy-flavour momentum correlation [Q, E3L although well defined in the 
massless limit, cannot be computed using the massless results of refs. |], ||, 

In this paper we describe a recently completed next-to-leading-order calculation 
of the heavy-flavour production cross section in e + e~ collisions, including quark mass 
effects. Very recently, two calculations have appeared in the literature that address 
the same problem || |7], §. They both use a slicing method in order to deal with 
infrared divergences. In our work, we preferred to use a subtraction method, since, in 
this way, we do not need to worry about taking the limit for small cutoff parameters^. 
We were able to perform a partial comparison of our result with that of ref . || , and 



found satisfactory agreement. In the older work of ref. [10], a calculation of the process 



e + e — > QQgg has been given, but virtual corrections to the process e + e — > QQg were 



not included. In ref. fll ], the NLO corrections to the production of a heavy quark 
pair plus a photon are given, including both real and virtual contributions. 

The paper is organized as follows. In Section [| we give a brief outline of the 
calculation. In Section ||| we introduce our kinematical definitions and conventions. In 
Section [| we present a somewhat detailed description of the calculation. In Section |5] 
we describe a few checks on our result. Finally, Section |6] contains some concluding 
remarks. 



Subtraction methods for the calculation of radiative corrections to e + e~ — * jets have been used in 
refs. p], 0], and they have also been successfully employed in the calculation of hadronic production 
processes. 
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2 Generalities 

We begin by showing in Fig. [I] the Feynman diagrams for a Born term (a), a virtual 
correction term (b) and two real next-to- leading contributions (c,d). Next-to- leading 




c d 

Figure 1 : Some of the diagrams contributing to the process Z/^—* QQ + X : 
a Born graph (a), a virtual graph (b), a real emission graph (c) and a real 
emission graph with light quarks in the final state ( d). 

corrections arise from the interference of the virtual graphs with the Born graphs, 
and from the square of the real graphs. Observe that we always deal with the cross 
section for the production of the heavy quark pair plus the emission of at least one 
extra particle (i.e. a gluon or a quark). The inclusion of virtual graphs with only a 
QQ pair in the final state is not needed if one computes three-jet-related quantities. 
Furthermore, since we deal with unoriented shape variables, the kinematics of these 
virtual graphs is fully specified, and in order to account for them it is enough to 
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include, in the final result, a two-body contribution normalized in such a way that 
one obtains the correct total heavy-flavour cross section at order a 2 s (see [12| and 
references therein). 

Virtual graphs, besides the usual ultraviolet divergences (which are removed by 
renormalization) , also have infrared and collinear divergences. These cancel when 
suitable infrared-safe final-state variables are considered. Our treatment of the in- 
frared cancellation is such that the final result is expressed as a partonic event gener- 
ator, in which pairs of weighted correlated events are produced. Shape variables are 
computed independently for each generated event, and histogrammed with the cor- 
responding weight. Infrared-safe shape variables give rise to finite distributions. No 
arbitrary cutoff is needed in this calculation in order to implement the cancellation 
of virtual and real infrared divergences, since this cancellation takes place between 
the two correlated events. Therefore, one does not have to worry about taking the 
limit for a vanishing soft cutoff. This method is similar to the one of ref. |J, which 
was used there to compute a large class of shape variable distributions for the LEP 
experiments. 

At next-to-leading order, several complications arise that must be considered. 
In fact, heavy flavours may also be produced by a gluon splitting mechanism, and 
diagrams with four heavy quarks in the final state are also present. Interferences 
between gluon splitting and direct production should also be considered. It is useful, 
however, to separate the various contributions in the following way. We examine 
each contribution in terms of cut Feynman graphs, which represent, individually, a 
single contribution to the cross section. We classify the cut graphs according to the 
following types: 

A) Contributions where the electroweak currents in the cut graphs are coupled to 
the same heavy-flavour loop, and there is a single QQ pair in the final state. 
These contributions are the most complex from the point of view of renormal- 
ization and soft and collinear divergences. They constitute the hard part of the 
calculation. They include graphs in which a pair of gluons or a pair of light 
quarks is present in the final state. We will call them A-type. We show some 
of them in Fig. ^|. 

B) Contributions where there are two QQ pairs in the final state. These include 
cut graphs with a single heavy-flavour loop coupled to the weak currents, graphs 
with two heavy-flavour loops, one of which is coupled to the weak currents, and 
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Figure 2: Some of the diagrams of A-type. Depending upon the cut, each 
graph represents a contribution coming from the square of the four-particle 
final state or from the interference between the tree-level graph with a virtual 
correction. 

graphs with two heavy-flavour loops, where each loop is coupled to one weak 
current. These contributions are finite, and their computation is a straightfor- 
ward algebraic problem. We will call them B-type. We have collected some of 
them in Fig. [|. 

C) Contributions where the electroweak currents are coupled to light quarks. Also 
these contributions are finite, and easy to compute. The heavy-flavour pair in 
the final state is generated by gluon splitting. We will call them C-type. 

D) Interference between terms in which the weak current is coupled to the heavy 
quarks and to quarks of different flavours. These terms have the structure of 
Fig. By Furry's theorem, they must vanish for vector currents. For axial 
currents, they cancel in pairs of up-type and down-type quarks, because they 
have opposite axial coupling. Thus, the up-quark contribution cancels with the 
down quark, and, if the charm mass is neglected, the charm contribution cancels 
with the strange. Only the graph with a top quark loop remains. We call these 
graphs D-type. 



E) Graphs with two heavy-flavour loop coupled to the weak current, one of which 
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/ w| -]- -f- jw 

I 

Figure 3: Some of the diagrams of B -type. 



is virtual. We call these graphs E-type. They pair naturally with the D-type 
graphs with the top loop, since in cases of practical interest the top loop is also 
virtual. 



Most of the following discussion will deal with A-type graphs, since the other cases 
are either straightforward, or they have already been considered in the literature. For 



example, graphs of B and C type have been computed in ref . |L0j , and graphs of type 
D and E have been considered in ref. Ul3|. There is, however, one extra contribution 



that should be considered together with the A-type graphs, that is to say, virtual 
graphs in which a heavy-flavour loop corrects the gluon propagator. These graphs 
are ultraviolet divergent, and so their inclusion is mandatory if one wants to have 
the complete cancellation of ultraviolet divergences after renormalization. We will 
discuss this contribution in detail when we deal with renormalization. 

A-type graphs contain ultraviolet, soft and collinear divergences that must be 
regulated. Soft divergences arise when, in addition to the basic QQg final state, an 
extra soft gluon is emitted, giving rise to a real soft divergent contribution. Collinear 
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Figure 4: Structure of cut graphs of D-type. 



singularities arise when the final-state gluon in the QQg process undergoes a real 
(virtual) splitting into either a gg pair or a qq pair. In the first case, the collinear 
gluon can also be soft, so that collinear singularities can overlap with soft singularities. 

In order to regularize all singularities we used dimensional regularization. At first 
sight, this procedure would seem in conflict with the presence of the axial coupling. 
In fact, for the class of graphs of A-type, there is a simple trick to avoid this problem. 
First of all, we notice that for unoriented shape variables the axial-vector interference 
cannot contribute. In fact, for the three-parton final state there are not enough 
momentum vectors to construct an invariant with an e symbol. For the four-parton 
final state one could in principle build such an invariant, but the cross section must 
be symmetric in the light parton momenta, so that such an invariant cannot survive. 
We then consider the case of a generic vector current coupled to two fermions with 
different masses mi and m 2 . One can then easily convince oneself that the case of 
the axial coupling can be obtained by setting mi = m and m 2 = — m, since one can 
turn — m into m by a chiral rotation. This procedure is bound to work if there are 
no anomalies involved in the calculation, and this is certainly the case for our A- 
type graphs. We will therefore proceed to compute the 0(a 2 s ) three- and four-body 
contributions in d = 4 — 2e dimensions. We will get a result of the form 




where the suffix 3 and 4 refers to the three- and four-body contributions. The ultra- 
violet, collinear and soft singularities will manifest themselves as single and double 
poles in 1/e in the three-body contribution, and as singularities arising from the 
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phase-space integration in the four-body contribution. We will call q the total in- 
coming invariant momentum and p, p' the momenta of the outgoing heavy quark and 
antiquark. Furthermore, we introduce the variables 

2q-p 2q-p' {q-p-p'f (n ^ 

xi = — J - , x 2 = — =— , y = 5 • ^ 2 - 3 ) 

gZ qZ gZ 

Here y characterizes the mass of the light system accompanying the heavy-quark pair. 
Thus, for Born and virtual graphs we will always have y — 0. 

In addition 

d<E>3 = dxi dx 2 Jz(xi, x 2 ) 

rf$4 = dx\ dx 2 dyd 2 Y J^xi, x 2 , y, Y) , (2.4) 

where Y represents the other two variables that are necessary to describe the four- 
body final state. In order to implement the cancellation of the soft and collinear 
singularities, we now imagine to compute some physical quantity G, function of the 
final-state variables. The reader may think of G as the combination of theta functions 
that characterize a histogram bin for some infrared-safe shape variable. In general the 
definition of G will be specified for any number of particles in the final state. Since we 
are only dealing with three- and four-parton final states, as far as we are concerned 
here, G is characterized by only two functions, G^(xi,x 2 ) and G^(x±,x 2 , y, Y). Soft 
and collinear finiteness of G will require that 

\im G^\x 1 ,x 2 ,y,Y) = G^\ Xl ,x 2 ) . (2.5) 
y — *o 

We will have 

da® 



Jda®G = J dx 1 dx 2 J 3 (x 1 ,x 2 )-^-G^(x 1 ,x 2 ) 

r da (2) 
+ J dx x dx 2 dy d 2 Y J 4 (x u x 2 , y, Y) G^( Xl , x 2 , y, Y) , (2.6) 

where each term on the right-hand side contains soft and collinear divergences that 
cancel in the sum. This formula will be rewritten in the following way: 

r r (da {2) r da {2) 1 

J da^G = j dx 1 dx 2 G^(x 1 ,x 2 )l^-J 3 (x 1 ,x 2 )+ j dyd 2 Y^-J A ( Xl ,x 2 ,y,Y)\ 

+ J dx 1 dx 2 dyd 2 YJ 4 (x u x 2 ,y,Y)l^-G^(x u x 2 ,y,Y)-^-G^(x 1 ,x 2 ) \ (2.7) 
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— (2) 

where a\ is chosen in such a way that it has the same soft and collinear singular 

(2) 

part as 04 , or, schematically 



lim^ = l. (2i 



v->° da 



rf$4 

The first term of eq. ( |2.7|) can be computed analitically. The 1/e single and double 
poles present in da^ /rf$ 3 all cancel with the poles arising from the dy d 2 Y integration 
of da^ /<i$4, and thus this term is finite. 

The second term in eq. ( |2.7|) , because of eqs. ( j2.5|) and ( plf) , has no soft and 



collinear singularities, and thus can be evaluated directly in four dimensions^. It 
is easy to see how the computation of this term can be implemented numerically. 
Assuming for simplicity that we can generate four-body configurations uniformly in 
the four-body phase space, to each four-body configuration x\,x 2 ,y,Y we associate 
two events: one four-body events with kinematics X\,x 2 ,y,Y and weig 
and one three-body event with kinematics x\, x 2 (and y = 0), and weight —da^/d^^. 
The computation of a shape variable using the above scheme reproduces exactly the 
second term of eq. (|2.7|). 

3 Kinematics 

3.1 Three-body kinematics 

We consider the following three-body process 

e+ (p'J + e~ (p e ) - Zh (q) - Q(p) + Q(p') + d(k) (3.1) 
where Q is the massive quark, and the momenta of the particles satisfy 

2 12 2 i 2 n 

p = p = m k = . 

Since we are interested in unoriented shape variables, we can express the three-body 
phase space in terms of two variables, which we choose to be 

2g • p 2q-p' 
xi = — 5—, x 2 = — 5—. (3.2) 
q z q z 



3 Observe that both eq. (2^) and eq. ( |2.8| ) must be satisfied in d dimensions in order for this 
argument to apply. 
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Defining in addition 



4m 2 

P = 

T 



(3.3) 



the three-body phase space in d = 4 — 2e dimensions takes the form 
(PS) (3) = H f 1 dx! r + dx 2 



4 (x\ - p) (x\ - p) - [x 2 g - {x\ - p) - [x\ - p) 



where: 
H = 

x 2 ± = 



16/T 



2c 



r(2 - 2e) 2 (4tt) 3 V q 2 
1 



4(l- Xl )+p 



2 (1 - Xl ) (2 - Xl )+ p {2- Xl ) ±2(1- Xl ) \jx\-p 



(3.4) 



(3.5) 



.(3.6) 



3.2 Four-body kinematics 

The four-body processes we are considering are 

e+ (p'J + e" (p e ) - Z/ 7 (g) - Q(p) + Q(p') + d{k) + <?(0 
e+ (p' e ) + e" (p e ) - Z/ 7 (?) - Q(p) + Q(p') + + g(Z) , 

where q is the massless quark. The momenta satisfy 

I 2 = k 2 = 

In the centre-of-mass system of the two massless particles, we have 

/ = l (1, . . . , sin0sin0, sin cos 0, cos 0) 
k = ko (1, . . • , — sin 9 sin 0, — sin 9 cos 0, — cos#) 



2/2 2 

p = p = m 



(3.7) 



P 



P 



Po i,...,o,o,Wi 



Po 



p'o 1,...,0, 



\ 



1 o sm 01 

Po 



1 o cos a 

Po 



where the dots indicate d — 3 equal and opposite components in the expression for / 
and k, and d — 3 zeros in the expression for p and p' . 
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To describe the unoriented four-body phase space, we need five independent vari- 
ables, which we choose to be 



x x = 
We thus have 
k = h = \/q 

and 



2q ■ p 



x 2 = 



2q-p' (fc + Q 2 

— — , y = — 2 — ' > ^- ( 3 - 8 ) 



2 



2 ' 



cos a = 



y(p-x 1 - x 2 ) + (1 - xi)(l - x 2 ) - 



Setting 
and defining 



yj(l -x 1 -yf -Apy yf(l -x 2 -yf - Apy 



v — -(1 — cos#) 



x fl — 


2- 


- Xi 


- x 2 




1 






y± = 


4 


±2^ 


M- 






2 


- p 


x 1+ = 


2- 


2- 





the four-particle phase space in d = A — 2e dimensions is 

.2 / 4?r \ * 



(p S )(4) 



5 



(47r) 2 r(l-e) \q 2 



-4 if 



X < 



X 



1 rx2+ ry+ rxi+ rX2- ry- 

dxi I dx 2 J dy + J dx\ \ dx 2 \ dy 



Vp Jx2 



y- 



rcgion I 



region II 



(3.10) 



(3.11) 



(3.12) 



{4 (x\ - p) (x\ - p) - [{x] - Ay) - {x\ -p)- (x\ - p)\y 

x j 1 dv[v{\ -v)]- e ^- r#(sin0)- 2e (3.13) 
J N$ Jo 



with 



N 4> = r d(f) (sin <f))- 2e = A e 7i 

J 



r(l-2e) 
P(l-e) 



(3.14) 
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region I 










region II 







Vp 



Figure 5: The two different areas in the %\-Xi plane correspond to the region 



I and to the region II of eq. (3. IS) 



As can be seen from Fig. [|, the integration region is split into two parts, one 
of which (region I) is characterized by the same X\ and x 2 integration limits as the 
three-body phase space. In this region, the variable y can reach 0, and therefore 
collinear and soft divergences arise. 

Sometimes we will need an analogous set of final-state variables, in which the 
role of p and p 1 are interchanged. The variable y remains the same, x\ and 22 are 
exchanged, and the other two variables, denoted by v' and (ft', are related to v and </> 
by the equations 



cos a — (sin 9 cos 6 sin a — 2v cos a 



1 — cos a — 2 (v — v' cos a) 



(3.15) 



COS( 



2 sin a Jv'(l — v') 

Exchanging instead the roles of I and k brings about the following transformations: 

V— >1— V, (j)^7T + (f), v' — > 1 — v' , (j)' — >7T + (j)'. (3.16) 
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4 Outline of the calculation 



The amplitude for the process can be written (up to an irrelevant phase) 



u(p e 



gl 



q 2 -M 2 + iTM z 



tvy* - 0.7V) (0| Jv(0)v Q ~ Jl(0)a Q \f) 



+ 9 2 ^f(cein(0\Jv(0)c Q \f) 



"(Pi) , 



(4.1) 



where / refers to states with four-momentum q. We use the notation 

= 9 

2 sin #w cos 0\y 

Vi = T 3i - 2ci sin 2 9 W 

where g is the electromagnetic coupling, T 3i is the third component of the (left) 
isospin of fermion i, Ci is its electric charge in units of the positron charge and #w is 
the Weinberg angle. Since we are interested in unoriented events, and following the 
assumptions described in Section ^ we can neglect the axial-vector interference in the 
square of the amplitude. From eq. (|4. 1|) we get the following cross section, averaged 
over the incoming electron beam direction 

N r 47Ttt 2 



da 



3q 2 



dT A 



p 2 {q 2 ) {v 2 + a 2 e )a 2 



+ dT v [p 2 {q 2 ) (vj + a 2 e ) v 2 Q + c 2 c 2 - 2 pi(g 2 ) v e v Q c e c Q 

where a is the electromagnetic coupling constant and N c = 3 is the number of colours. 
In addition 

P2(q 2 ) 

We have also defined 



q 2 (m| - q 2 ) 



4 sin 2 9 W cos 2 6» w (m| - g 2 ) 2 + m% T% ' 



4sin 2 6> w cos 2 6'w/ (jn% - q 2 ) 2 + m% T% 



M 



V/A 



(/») 

V/A 



EM 

n 

2tt 



V/A "rn 



+ q 2 



(0\J v/A (0)\f n )(f n \J v/A (0)\0) 



(4.2) 



-14- 



where d<j) n represents the n-body phase space, and f n represents an n-body final 
state. The q^q v term in the projector in eq. ( |4.2| ) is, of course, irrelevant for the 
vector current component, but it should be kept for the axial current when the quark 
mass is non-zero. 

In the following we will be interested in strong corrections up to the second order, 
and into the final states: QQ, QQg, QQgg and QQqq. We will use the following 
simplified notation: 

(2) — 

• M. V , A for the QQ Born term 

• M-v)a or to indicate the three-body QQg, 0(a s ) term 

• M-v)a or M v to indicate the three-body QQg, 0{a 2 s ) term 

• Myf A or M gg for the four-body QQgg, 0(a 2 s ) term 

• Myj\ or Mgg for the four-body QQqq, 0(a 2 s ) term, 

and equivalent ones for the dT v /A terms. 

We will drop the V/A suffix when not referring specifically to the axial or vector 
contribution. 



4.1 QQ cross section 

From the amplitude 

u(p)T v/A v(p') , (4.3) 

where T v = 7 M and T A = 7 M 7 5 , we obtain the two-body cross section at zeroth order 
in a s . We get 



where p is defined in ( |3.3| ) and 

P = y/T^~p. (4.5) 
Multiplying eq. (fOj) by the 2-body phase space (3/(8%) we get the zeroth-order total 



cross section 



7f=/3(l + 0, Tf=/? 3 . (4.6) 



Thus, our choice for the normalization factor in eq. ( (4.2|) is such that, in the massless 



limit, T v 2) = T A 2) = 1. 
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4.2 QQg cross section at order a s 



This is obtained starting from the amplitude 

4 + U + m . „ u 



We define 



u(p) 



1 (p + k) 2 - m 2 V ' A V ' A (p - g) 2 - m 2 r 



^ q 2 J Z^^V/A^V/A 



v(p) ■ (4.7) 



(4.8) 



where the sum refers to the spin of the fermions in the final state. 
The sum over the gluon polarization is 



My i A = -g aa , M 



era 

V/A • 



(4.9) 



We will need M v / A in d = 4 — 2e dimensions. We have 



M 



1 T" Q 



V 



+ 



16 



(l-xi)(l-x 2 ) (l-X!) 2 (l-a; 2 ) 2 \g 2 
-3 (a; 2 + xj) - 8 (1 - Xl )(l - x 2 ) + 2 
16e 



2xix 2 (xi + x 2 ) 
32 (m 2 



(I - Xl ) 2 (l - x 2 ) 2 \q 2 



x, 



(1 - Xi)(l - x 2 ) 



x x + x 2 + (1 - Xi)(l - x 2 ) + x g — 1 



m 



.17 



g 2 I (1 - a; 1 )(l - x 2 ) 



+ 



xr 



{l-x 1 ){l-x 2 ) 9 



(4.10) 



and 



Mi 



T2 I ,-y> 2 
1 ~~ 2 



+ 



16 



m 



(1 - m)(l - x 2 ) (1 - xi) 2 (l - x 2 ) 2 V q 2 
-llxix 2 (xi + x 2 ) + 8(x 2 + x 2 ) + x\{x 2 — 1) + x 2 (xi — 1) 



-12(xi + x 2 - 2xix 2 ) 



+2x(x 2 2 + 4 (1 - e) + 



64 



m 



2\ 2 



(l-Xi) 2 (l-x 2 ) 2 \ g 2 



(1 - e) 



(1 - Xi)(l - x 2 ) 



X 



9 ■ 



(4.11) 



where X\,x 2 and x g are defined by ( |3.8| ) and ( ft.llp . The three-body, order-a s cross 
section is given by 



Mty A = ^C F g 2 s ^M v/A 
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where \i is the mass parameter of dimensional regularization and Cf = * = I ^ or 
an SU(3) gauge theory. 

We introduce a unit, purely space-like vector j lying in the event plane (i.e. the 
plane defined by p, p' and k) , and orthogonal to k 

j.q = 0, j-k = 0, f = -l. (4.12) 

Myj' A has the general form 

Myj' A = My- /A g aa> + M v/A ff' + terms involving q or k . (4.13) 

In the following, we will need My, A , but only in four dimensions 

4xix 2 (xi + x 2 ) - p(xi + x 2 ) 2 - A(x\ + x\) 



Ml 2CV/A 



V/A (l-xi) 2 (l-x 2 ) 2 

-12x1X2 + 4(p + 2)(xi + x 2 ) - 4(p + 1) 

where 



(4.14) 



c v =p + 2 c A = -2(p-l). (4.15) 

We also define, consistently with our previous notation 

M%? = ^ C F fifPM# A , M^ /j = ^ C F gl^M^ A . (4.16) 



In the cases when the V/A suffix needs not be specified, we will simply write Ai™' , 
Mi and M{. 



4.3 Virtual contributions 

Corrections to the three-jet decay rate to order a 2 s come from the interference of 
the one-loop graphs with the tree-level ones. These terms have been computed in 
d = 4 — 2e dimensions. The algebra has been carried out in a straightforward way, 
using a MACSYMA program, which reduces the original Feynman graphs to a linear 
combination of scalar, one-loop integrals. The scalar integrals have been computed 
analytically. Their values are listed in Appendix |Fj. Loop corrections to on-shell 
external lines require particular attention. First of all, gluon and light fermions self- 
energy corrections to external gluon lines vanish in dimensional regularization. Only 
the corrections coming from a heavy-flavour loop need be considered. We proceed 
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as follows. We compute the self-energy correction for a gluon propagator of small 
virtuality. We obtain, for the corrected propagator, 

- N,T F £ ( 41 ' i_ -'fa--MW t2 > + ( ^p-, (4. 17) 



k 2 s \m 2 J 3e k 2 k 2 

where 

N e = ^T{l + e) (4.18) 

and the colour factor Tp = 1/2. 

From this equation we immediately infer that the contribution to A4 V coming from 
the self-energy correction to the external gluon line amounts to 



' W lS^ XM - (419) 



A similar consideration applies to the self-energy corrections to heavy-flavour external 
lines. In this case one finds 

2/VW3 A i i .. i 0{p 2 -m 2 ) , Ann . 
N e C F g 2 s - + 4 - - i5m- + J - , (4.20) 



]/> — m s \m 2 ) Ve J $ — m $ — m $ — m ]/> — m 

where 

5m = N e C F g 2 s (j^ m + 4) . (4.21) 

The infinite mass correction should be removed by the mass counterterm. We define 
the Feynman rule for the mass counterterm to be given by an insertion of —im c in 
the fermion propagator, where 



m r = —5m 



-A e C^^Vm0+4) . (4.22) 



This precisely cancels the 5m term in eq. (|4.20| ), so that the pole of the propagator is 
not displaced by radiative corrections, and m corresponds to the pole mass definition. 
Thus, the effect of the fermion self-energy correction to an external line, including 
the effect of the mass counterterm, is given by 

- N < C '£(£)\- e +i ) XMb - (423) 

To complete the computation of virtual corrections, the diagrams with a mass coun- 
terterm insertion in internal fermion lines should also be included. After that, charge 
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renormalization is all that is needed, since we are computing a physical cross section. 



We carry out the charge renormalization in the mixed scheme of ref. [I4|, in which 
the light flavours are subtracted in the MS scheme, while the heavy-flavour loops 
are subtracted at zero momentum. In this scheme the heavy flavour decouples at low 
energy. The prescription for charge renormalization in this scheme is 



a. 



a. 



l + g 2 s N € 



3e 



T F n\ { 



11 

37 



C A 



JL 

m? 



3e~ 



(4.24) 



where Ca = N c = 3 for an SU(3) gauge theory. It amounts to adding the following 
correction to our virtual term 

An, 



UK) It, 



3e" 



x Mi 



(4.25) 



Observe that in this scheme the term corresponding to the heavy-flavour loop com- 
pensates exactly the self-energy correction to the external gluon line, coming from 
the heavy-flavour loop. This is easily understood: the final-state gluon is on the mass 
shell, so it is effectively renormalized at zero momentum by the heavy quark loop, 
and thus decoupling applies. We can now resume the combined effect of external line 
corrections and renormalization to be included with M v 




2C F - + 4 




x Mi 



(4.26) 



The factor of 2 in front of the fermion external line corrections is to account for the 
two fermion lines. 



4.4 Soft and collinear limit of the QQgg and QQqq 
cross sections 

Here we derive an expression for the singular part of the four-body cross section, 
valid in both the collinear and the soft limit. These limits are both characterized by 
y — > 0, except that in the soft limit, at the same time v — > (I — > 0) or v — ► 1 (k — > 0). 
We will focus our discussion on the QQgg final state. The other process QQqq is much 
simpler, since only collinear singularities are present there. Since the same formulae 
apply irrespective of the vector or axial case, we will always drop the V/A suffix. 

We begin with the soft singularities of M gg . They are given by eq. ( |D.5| ), which 
we now rewrite 

+ 2 (r ° A \ p ' p ' 

\ F 2 J (p ■ l) (pf ■ i) 



Mlf ~ g\^ \C A 
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p ■ k p' ■ k 



(p-l)(k-l) (p'-l)(k-l) 



-19- 



m m 
+ 



+ (k <-> I) } x M 



b ■ 



(4.27) 



(p-l) 2 (p'-l) 2 _ 

From Section we can derive an approximation of the scalar products in the limit 
of I soft 

p-k 2h 1 



(p-l)(k-l) q 2 y[y + hv] 



E pM i(x 1 ,x 2 ,y,v) 



p ■ p 



K 1 



(p ■ l)(p' ■ I) m 2 y + hv y — c cos 4>y/y~y/v + g ■ 



E p ,p'-i(x 1 ,x 2 ,y,v, 



m 2 Ah 



(p ■ I) 2 q 2 [y + hv} 2 
where 

,2 



E PiP; i(x 1 ,x 2 ,y,v) , 



T 



(1 - x x ){\ - x 2 ) 
2 

(1 -xi)(l -x 2 ) 
1 - x 2 4 



Xi + a; 2 - 1 5- 



l — x 2 1 — Xi 



+ 



2 r 



1 — Xi 1 — 2 2 



1 — Xi 



1 - x 2 



1 — xi 6 
2v^ 



Xi + x 2 - 1 



2m 2 



*=6" 

We will also need analogous formulae in the variables in which the roles of p and p' 
are interchanged. We have 

p' ■ k 



E 'p',k;l( x ^ X 2, V, v ') = Ep, k -i(x 2 , Xi, y, v') 



p ■ p 



(p-l)(p'-l) 



E P ,p>;i( x i, x 2, y> v 'i = E p,p'A x 2, xi, y, v', (f)') 



m 



(p> ■ iy 



E 'p',p'-i( x ^ X2 ' V> v> ) = e p',p'A x 2, xi,y, v') . 
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Soft factors for the case when k is soft are instead obtained from the above using 
eqs. ( p. 16 ). For example 



E Pjl . k (x 1 ,x 2 ,y,v) = E Ptk .i(x 1 ,x 2 ,y, 1 - v) , 
E P rf ;k (xi,X2, y, v) = E P:P ,. l (x 1 , x 2 , y, 1 - v, <f> + n) . 

We can now write down our approximate soft cross section. We have 



M 



soft 
99 



9s ^ 2e jc 'a E pM i + E' p , Ml + E p j. k + E' p , Mh 



+(Cf — Ca/2) E PtP i;i + E pp , ;l + E PiP /, k + E' pp ,. h 



Ep, p -i + E p , tpl;l + E P:P . k + E' p ,p,. h 



x Mi 



(4.28) 



(4.29) 



The soft cross section written in this way is symmetric under the interchange of k 
and Z, and of p and p' . 

The collinear part of the cross section receives contributions from both the gg and 
the qq final state. For the gg contribution, according to eq. (P.10| ), we can write the 
collinear part 



2 2e^A 

g s v -5- 

q 2 y 



1 1 

-2 + - + - + z(l-z) 



z 1 — z 



9o 



2z(l-z)(l-e) 



k±ak ±(T ' g± c 



k* 



2 - 2e 



x .M 



6 5 



(4.30) 



where z is the momentum fraction of I versus I + k in the collinear limit. It can be 
chosen to be equal to v or to v'. 

The perpendicular direction refers instead to a direction orthogonal to I + k in 
the centre-of-mass system and in the collinear limit. Using eq. ( |4.13|) , the azimuth- 
dependent term of ( |4.30| ) becomes 



g 2 s ^^l-2z(l-z)(l-e) 



q 2 y 



Mt + Ml—p 



(k±-j) 2 Mi{2-2e)-M\ 



2 - 2e 



2 2<e 

9 s V 



4Ca 
q 2 y 



- z(l - z)M{ 



^2(1 - e ) + 1 



(4.31) 



It is now easy to show that, in the collinear limit, (k± ■ j) /kj_—> — cos 
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Part of the collinear singularities are already contained in the soft-limit expression. 
In fact, for y — > at v fixed, we have 



E 



p,k;l 



q 2 yv 



E. 



p,l;k 



q 2 y(l - v) 



(4.32) 



Thus, the 1/z and 1/(1 — z) terms in the collinear limit formula (}4.30|) should not 
be included, since they are already present in the soft term. We thus arrive at the 
following expression for the collinear term to be added to the soft term 



99 



Mi 



9 S V ^^{M b 



wcoii .:i.;i< i 



q 2 y 



v(l - v) + v'(l - v') 



1 - v) (2(1 - e) cos 2 (j) - l) + v'(l - v') (2(1 - e) cos 2 <f)' - l)] 1 (4.33) 



An analogous procedure yields an expression for the collinear part of M q q (see 
eq. (ED) 



1 



-qq 



q 2 y 



Ml r 



v(l — v) ( 2 COS — 



4(1 -e) 
1 



v ,2 + (l-v'Y + v 2 + (l-vy-2e 

1 



1 - e 



+ v'(l-v') 2cos : 



-,2 i/ 



1 - e 



(4.34) 



The expressions M°°g l , M^f 1 and -Mg° ft depend upon x\ and a; 2 via Mb and .M^. 
These expressions are meaningful only if X\ and x 2 belong to the domain of the 
three-body phase space. We thus define 



M 



mi 



M q q * 



M s g f + M c g f) 9 3 (x 1 ,x 2 ), 
AC- 11 9 3 (x u x 2 ), 



(4.35) 



where the 9% function is precisely defined to be zero when x± and x 2 are outside the 
three-body phase-space region. More specifically, using the definitions of eqs. (|3.12| ) 

9 3 (x 1 , x 2 ) = 9(1 - xi) 9(xi - y/p) 9(x 2+ - x 2 ) 9(x 2 - x 2 J) . (4.36) 

We are now in a position to specify the subtraction procedure outlined in Section |[ 
Our expression for the second-order contribution to an infrared- and collinear-safe 
quantity G is given by 



- / M gg (x 1 ,x 2 ,y,v,(f))G(x 1 ,x 2 ,y } v } (j)) u <>i 



+ J M qq -(x 1 ,x 2 ,y, v,<j>) G(x 1 ,x 2 ,y,v,(p) dfa + J M v (x 1 ,x 2 )G(x i ,x 2 ) 
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where all quantities are computed in d = A — 2e dimensions. The factor 1/2 in front 
of the gg contribution accounts for the two identical gluons in the final state. We 
rewrite the above expression as 
1 



M gg (xi,x 2 , y, v, <t>)G(xi, x 2 , y, v, 4>) - M gg (x l ,x 2 , y, v, 4>)G{x x , x 2 )j d$ 4 

[Mqq(x 1 ,x 2 , y, v, 4>)G{x u x 2 , y, v, 0) - M qq {xi,x 2 , y, v, 4>)G(x h x 2 )j d$ 4 

+ J (M v {x u x 2 ) + Mi{xx,x 2 )) G{xi,x 2 ) d$ 3 (4.37) 
where we have defined 

M i {x 1 , x 2 ) = ^ J Mggixx, x 2 , y, v, <p)d® 4/3 + J M qq -{x 1 , x 2 , y, v, 0)d$ 4/3 (4.38) 
and rf$4/ 3 is defined by 

d$ 4 9 3 (x 1 ,x 2 ) = d$ 4/3 d$ 3 . (4.39) 

An explicit expression for g?$ 4 / 3 can be obtained from eqs. ([3.13|) and 03.4|) . We first 
notice that the four-body phase space is almost proportional to the three-body phase 
space, except for the ratio 

(A(x\- p)(xl- p) - \(x 2 - Ay) - (x 2 - p)-(x 2 2 - p)] 2 \~ e s 

2 W 2 r r 2 , 2 \ , 2 M2 = 1 + ° f £ ' 4 " 40 

V A{x\- p){x\- p)-[x\-{x\- p)-[x\- P)} 2 J 
On the other hand, terms of order ye can be neglected, since they cannot generate 
infrared singularities, because of the y factor, and therefore they can only produce 
terms of order e. Thus we can write 

d<S> 4/3 = N e R e q 2 q- 2 * f + dyy'* f dv [v(l - v)}~ e -L f # (sin 0)~ 2e (4.41) 
Jo Jo JMtp Jo 

or the analogous one in the v ', <$' variables. The normalization factor N e is defined 
in Q4TT8D , while 

Re = -, ^ r = 1 - — + 0(e 3 ) . 

r(l + e)r(l-e) 6 V ; 

Since we are free to choose the set of variables we prefer in the <i$4/3 integration, it 
is easy to see that the J\4i(xi,x 2 ) term reduces to 

Mi(xi,x 2 ) = g y* [\l^(v(l-v)-2) 

J 12 q 2 y 

V + (l-t;) 2 - 



2nuTp 1 

+ 



q 2 y 1 - e 



+ \[AC A E pMl + A{C F - C A /2)E PjP , l - AC F E P ^} 1 x M b rf$ 4/3 
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where the term proportional to }A? has been dropped, since it vanishes in d = 4 — 2e 
dimensions, after the azimuthal integration. 

Furthermore, the remaining collinear term is easily integrated. We define 
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V+ _ /•! 1 

dyy e / dv[v(l-v)\ € — 
n Jo iV^ Jo 



smi 



\-2e 



[«(!-«) -2] 



1 



[1 - elog(y+)] 



11 _ 67 
If ~ 18* 



+ 0(t 



(4.42) 



and 



rcoll 
l qq 



/ A/y- / *Ki-«rv / < 

jo Jo iV^ Jo 

~[l-clog(y + )] (| + ye)+0(e) 



sini 



-2c 



lt; 2 + (l-t)) 2 -e 
y 1 — e 



For the integrals of the soft term we define 



I p ,k;i = q 2 / e / dv[v(l-v)) 



1 

/o Jo " Nfj, Jo 

and the analogous ones for J PiP; ; and In this way 



sm< 



-2e 77. 

&p,k;l 



(4.43) 



^p,fc ; « — 2/i 7i J PiP; i — 4/i I 2 1 P,p';l 

where the values of Ii, J 2 and I 3 are collected in Appendix |F|. 
Our final expression for j\4i(xi,X2) is therefore 

Q 



I-n -n' d — K I3 



M i {x 1 ,x 2 ) 



N e R t g 2 s (^) \ 2C A I c g f + 2 n H T F I c g f 



+ 



2CAl P ,k-,i + 2(Cf — Ca/2)I p , p ';1 — 2CfI p , P ;i 



x M b . 



5 Checks of the calculation 



Several checks have been performed to control the correctness of our results. 

1. The divergences coming from UV and IR poles all cancel. 

2. The full calculation, asm^O, agrees with the massless result of ref. [|TJ. 
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3. Our four- dimensional matrix elements for the processes e + e — > Z/7 — > QQgg 



and e + e~ — > Z/7 — > QQQQ agree with ref. [10]. Furthermore, the soft and 
collinear limits of the four-body matrix elements for the process Zj^^QQ 
plus two light partons are correctly given by formulae fl4.35|) . 



4. Near the production threshold, we should recover the Coulomb singularity. If (3 
is the velocity of the two massive quarks in the fermion centre-of-mass system, 
then (see ref. (ffBl) 



da^J A (x h x 2 ) ^ (c F ~ ^) da i y ) /A (x 1 ,x 2 ) . (5.1) 

By evaluating {p+p') 2 in the centre of mass of the two massive quarks, for small 
(3, we get 

r / ™ \i2 

= (q-k) 2 = q 2 (x 1 +x 2 - 1) • (5.2) 



2[m + -^ + 0{n 



{p + p') 2 

Choosing for example x\ = x 2 we have 

-(l+p + p/3 



X\ = x 2 



By letting j3 get smaller and smaller we have checked that the behaviour of the 
virtual differential cross section is in agreement with eq. (|5.1|). 



A further check is described in detail in ref. 



6 Conclusion 

In this paper we have described a next-to-leading-order calculation of the heavy- 
flavour production cross section in e + e~ collisions, including quark mass effects. Some 
applications of our calculation have appeared in the literature 0, [0. 

We have used a subtraction method instead of a slicing method, in order to avoid 
having to worry about taking the limit of some small parameters. We have performed 
several checks on the correctness of our results. Among them, the small mass limit 
of the energy-energy correlation is of particular significance, since, for this quantity, 
some discrepancies among different approaches are still present (see ref. [I7|). 
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Appendix A: Phase space for four massive quarks 

For completeness, we describe in this appendix the phase space for four massive 
quarks in the final state. The process is 

e+ (p' e ) + e~ (p e ) -> Zh (q) -> Q(p) + Q(p') + Q(r) + W) , (A.1) 

where 



r 2 _ r /2 



2/2 2 

p = p = m 



The four-body phase space is obtained with a procedure similar to the one given in 
Section ^|, with the simplification that now the entire cross section has no infrared 
or collinear divergences, so that we can put ourselves directly in d = 4 dimensions 
and we do not need to divide the phase-space region into two different pieces. In the 
centre-of-mass frame of one heavy quark-antiquark pair we have 

r = (ro, |r| sin 9 sin 0, |r| sin 9 cos 0, |r| cos#) 



(r , — |r| sin 9 sin 0, — |r| sin^cos^, — |r| cost 



p = p I 1,0,0, y i - p 2 



P'=P' I 1-U. 
where 



a 1 72 sina, 

\ Po 



, 1 cos a 

\ Po 



y 



(r + r') 



A2 



fo = yjq A 



q« ' " 2 

and p ; Po an( l COSQ; are given by ( ^.9| ) and ( |3.10| ), while |r| = y — m- 
The four-body phase space is given by 



(PS) 



where 



(4)__9 



2/H 



(1 , dyJl / rfxi / cfx 2 dv 

(Ait) Jp \ y J^p Js 2 _ Jo Jo 



Xi+ = 1 - 2/ - v 7 ^ 
1 



^2± 



4(1 -xi) + p 



(2-x 1 )(2 + p-2y-2x 1 ) 



(A.2) 



±2y/(xl-p) [{xr-l+yf-py} 



(A.3) 



-26- 



A statistical factor l/(2!2!) = 1/4 must be supplied to (1A.2Q , because of the presence 
of two pairs of identical particles in the final state. 



Appendix B: Collinear limit for g^gg splitting 



a 




v 



Figure 6: Gluon splitting 



In this appendix we will derive the singular part of the square of the invariant am- 
plitude when two collinear gluons are produced. In the collinear limit, the amplitude 
for the emission of two gluons in the final state can be decomposed into two parts: 
the first one contains the graphs where the two gluons are emitted by a single virtual 
one (see Fig. ||), and the other one contains all the other graphs 

A ah = {aI{1 + k) + ° (-g s ) f^V^i-k, -I, k + l)+ K%} e,(k) c„(J) , 

(B.l) 

where a and b are the colour indices of the final gluons, P is the spin projector of the 
gluon propagator, g s is the strong coupling constant, f abc are the structure constants 
of the SU(3) gauge group, e and e are the polarization vectors of the final gluons, and 
p^7 j g L oren tz part of the three-gluon vertex 

T^{-k, -l, k + l) = {-k + l)i gi" + {-21 - ky g v i + (2k + l) v g^ . (B.2) 

Only the first term of ( p.l| ) is singular in the collinear limit. We want to stress the 
fact that this term is singular in the soft limit too. Therefore one has to be careful, 
when considering the soft and collinear limit of the square amplitude, not to include 
this contribution twice. 
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We introduce two light-like vectors 



7/ 



C X 



k + l , k + l 

1 k + l 



k + l 



k + l 



(B-3) 



and choose c = 1/4, so that 2t ■ r] = 1. We then decompose 

Z" + Jfe' 4 = # i + f jj" , 

where 



(B.4) 



e = (J + *) 2 = ?V 

We will work in the light-cone gauge, characterized by the light-like vector rj, because, 
in this gauge (as we will see), the interference of the divergent term of ( |B.1[ ) and of 
the finite term 1Z does not contribute to the singular part. 

The gluon spin projector is then written 

if]) 1 + rpp c 



P a \p) = -g al + 



7] ■ p 



(B.5) 



We write I and k as 



(B.6) 



k» = vt>* + gr) tl + k l t 

p = (l - v ) t" + £V - k± , 

with k± such that t ■ k± = rj ■ k± = 0. By imposing that k 2 = I 2 = and that 
(I + k) 2 = q 2 y we have 



ki 



-v(l-v)q 2 y £ = {l- v )q 2 y £' = vq 2 y 



From ( |B.4| ) and ( |B.6| ) we finally obtain 

1 



k" 



1-v 



vl fl + (l-2v)q 2 yr]^ + k^ 



1 



(B.7) 



(1 - v) k>* - (1 - 2 w)g 2 ?/^ - fc£ 



Considering that 



(fc + 7 ^ 7 = 2q 2 yr]a k^{k) = /^(Z) = , 
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with the help of eq. (|B.7|) we can rewrite the amplitude ( ]B.1|) 

jab 



X 



q 2 y 

-2 kl gT + ~k± g ul + — ^ — kl + 
f 1 — v 



+ n%\e»{k)e v {l). 



Observe that the first term is of order , so that a singularity with strength 1/y 

can arise only from the square of the first term, and the interference term does not 
contribute. Furthermore, we can now substitute 

Alii + k)->A a c (t) = tree-level amplitude 

ifF + rpt c 



9 U, + 



7] ■ t 



-gT 



Remembering that r\ ■ e = r\ ■ e = 0, we can write A in the form 



92 



A ab = A ca (t)^ftf 
TV 



abc 



v 1 



e,(k)e u (l), (B.8) 



where only the term contributing to the singularity has been kept. 

By squaring the amplitude and summing over the colours and spins of the final 
gluons, we obtain, for the collinear singular part 



gg g2 y 



-2 + I + J- 
v 1 — v 



2v{l-v)(l-e) 



+ v(l-v) g 
k± a k± a > g Laa 



A°(t)A*/(t) 



(B.9) 



k\ 2 - 2e 

where we have used the gauge invariance t a A"{t) = to write the following identity 

A:(t)Af(t)g ±aa ,=A:(t)Af(t)g aa ,. 



The first term of ( |B.9| ) is recognized to be the Altarelli-Parisi splitting function for 
the gluon-gluon process, in d = 4 — 2e dimensions. The second term vanishes after 
azimuthal average in 4 — 2e dimensions. 

Coming now to our problem, we can further specify the structure of A a c {t)A*f (t). 
In fact, by using eq. ( |4.16|) , we can write ( P-9|) in the following form 



A ^col 2 2e 

M m = g s n 



4CU 
q 2 y 



-2 + - + - 

v 1 



2v(l-v)(l-e) 



- v 
k± a k 



9o 



+ v (1 — V 

g^crcr 1 

~~ 2 - 2e 



x M 



era 
b 



(B.10) 



-29- 



Appendix C: Collinear limit for g^qq splitting 



Following the same steps as in the previous appendix, we can give the approxi- 
mation of the square of the amplitude in the limit of a collinear couple of massless 
quark-antiquark. The invariant amplitude is 

A = A«(k + l) iP ^ k + l)6C "- 



u(k) (-tg s yt c ')v(l) 



(k + l) 2 

where P is given by (|B.5|) and t c are the generators of SU(3) gauge symmetry. By 
squaring this amplitude and summing over the spins and colours of the final quarks, 
we obtain 



M 



col 

qq 



T F gl 
<fy 2 



A a c {k + l)A?'(k + l)P pa (k + l)P p/(Tl {k + 1) Tr (liYlY') , 



where t c are normalized such that Tr (t a t b ) = T F 5 ab . 

Considering now eqs. ( |B.6[ ), we see that, in the collinear limit, the trace is of the 
order of y, so that the singular part can be obtained by putting y = in the rest of 
the numerator 

T F g 2 s 



M 



col 

qq 



■Al{t)Af{t) g^g^e' Tr {Hi P h P ' 



where we have used the definition of t given in eq. ([B.4[) . Evaluating the trace and 
keeping in the numerator only the terms proportional to y, we obtain 

Tf9 2 s 
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■qq 

that is 
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2 - 2e 
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v 2 + (l-v) 2 -e 
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(C.l) 



k\ 2 - 2e 

Here again we can recognize the Altarelli-Parisi kernel for g^qq splitting. 

As done before for eq. ( p.9| ) , we can specify this formula to the problem we are 
studying. With the same substitutions made to go from eq. ( |B.9|) to eq. (p.lOj ), we 
can write 

4 Tp f 1 



M% = gl^ 



q 2 y 



2 - 2e 
+ 2v(l- 



+ (1- 

k± a k A 

k 2 



9o 
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(C.2) 
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Appendix D: Soft limit for the invariant amplitude QQgg 

In this appendix we will derive the divergent part of the invariant amplitude for 
the process 

Z/l(q) - Q( P ) + W) + 9(k) + 9(1) (D.l) 

in the limit when the momentum I of the gluon is soft. A soft singularity appears 
only if the soft gluon is emitted from one of the external legs. If the emitting external 
particle is the gluon, the amplitude of the process, in the Feynman gauge, is 

A% {s) = A%(1 + k) ^jy 2 (-g s ) f abc F^(-k, -I, k + I) e,(k) e v {l) , 

where we have added to eq. ( |B.1| ) the colour indices i,j of the produced quarks. 
As I goes to zero, this term develops a singularity. By using the gauge condition 
k a A l J a (k) = and the transversality k^e^ik) = 0, we can write the amplitude as 

A-- is) = g s f abc ^— A%(k) e a (k) e„(Z) + non-singular terms. (D.2) 

Similarly, if we consider the emission of a soft gluon of colour index b from an 
external quark leg with momentum p and colour index i, that is 

Qn(p + l)^Qi(p)+g b (l), 

we can write the invariant amplitude 

Af Q) = u{p){-ig sl v t\ n ) *_ m AT,(P + I) e„(*) e v {l) , 

where A refers to the rest of the process from which the quark external line takes 
origin. 

In the limit of I going to zero, we can rewrite this amplitude as 



b P 

Aij = g s — 7 1\ A^ip) ejk) eJl) + non-singular terms, (D.3) 
p ■ I " 



where we have defined A^{p) = u(p)A^'j(p)- 

In the same way, we can obtain the limit of the amplitude for the soft emission 
from an antiquark with momentum p' and colour index j 

^afe(Q) = _^ __ AZ{p) t b nj e^k) e v {l) + non-singular terms. (D.4) 
p ■ l 
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Considering that A^J = t^A a } where A a does not contain any colour element, and 
the similar ones for eqs. ( D.3 ) and ( OA ), we can sum the three amplitudes to obtain 



-A-ij — 9s ^if k . p . I p' ■ I | ^ e v{f) 

where we have disregarded the non-singular terms. 

By squaring the amplitude and summing over the spins and colours of the final 
gluons and quarks, we have 

p ■ k p' ■ k 



M s ^(l)=gi^{-C 



+ 



(p-l)(k-l) ( P '-l)(k-l) 

,2 



+ 



P • p 



m 



+ 



m 



(p-l) 2 (p'-l) 2 



x Ml 



(D.5) 



2 J (p-l)(p'-l) 
where we have made use of eq. ( |4.16|) . 

The same result applies in the case of k soft, once the interchange I <-> k is made. 



Appendix E: One-loop scalar integrals 

We can classify the different types of scalar integrals according to the number of 
massive propagators in the loop and according to the "shape" of the loop: boxes (B) 
and triangles (T). We introduce the following kinematical invariants 



01 = (Q — P) 2 - m 2 = q 2 (l - x 2 ) 

02 = (q - P) 2 -m 2 = q 2 (l - x x ) 

03 = (q ~ k) 2 = q 2 (l - x g ) , 
where X\ and x 2 are defined by ( |3.8|) and x g by (|3.11|) , and 



(E.l) 



P± = 2 



a.\ ± \ la 2 — 



Am 2 



with : ai = 1 \ 

q 2 
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r/ ± = i(A'±A). 

Here we also give the absorptive parts of the integrals, although they do not contribute 
to the cross section. The integrals are computed in d = 4 — 2e dimensions. Terms of 
order e or higher have been dropped. 

Defining the dilogarithm function as 

log(l - z) 



Li 2 (x) = - f 
Jo 

and collecting the same factor 



dz 



N(e) 



16 7T 2 

in front of each expression, we obtain 
d d l 1 



(4tt) £ r(l + e) =iN € 



B- 



2m 



B 



3m 



rpq 

1 2m 



r dH 1 1 1 1 

= J (2n) d I 2 (I - k) 2 (l + p-q) 2 -m 2 (l+p) 2 -m 2 

, r . . / n\-e 1 f 1 1 /, m 2 , m 2 \ n , m 2 , m 2 

= iV(e) m + - log— + log— +21og— log — 

v ' o x o 2 [ e e \ o"i 02/ o"i ^2 



- ^7r 2 - log 2 ^± + 2m 
d d Z 1 1 



1 2 2 

1 to to rt + 

- + log h log h log — 

e 0\ o 2 A 



(E.2) 



1 



1 



(27r) d I 2 (I + p) 2 — m 2 (I + p — g) 2 — m 2 (7— p') 2 — m 2 



= iV(e) (to 



2\-e 



1 



1, e 



m 



o 2 o 3 A' I e 



log 7- + 2 log — + lo. 



0"2 



to 



&3 



V+ 



- 2 Li 2 (£_) - log 2 £_ - 2 Li 2 I -—J - 2 Li 2 ^ 



log : 



1 2 ^+ 

log — 

V+ 



+ 2Li 2 



, 7 K 



+2Li2 (_^ +log2 (_^ +2Li2 (_%)_ 



7T 



?7_ 



+ 2 log ( ] log ^ + log 2 e+ - 2 log A' log£_ + log 2 A' 



d d l 1 



V- 



1 a 2 £ 

- + 2 log — + 4 log r] + - 2 log A' + 2 log 

e 02 ^ 



(E.3) 



/ 



1 



1 



(2n) d I 2 (l-p') 2 -m 2 (l+p + k) 2 -m 2 
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N(e) (m 2 ) 



2\-e 



+ Li 2 



Lis 



+ Li, 



' K - p+ 

^+ -p- 
, i-p- 

>_-A_ 



^ Li2 fi_r 

2 -^#l V p+, 



Li, 



P+ 



A+ -p+ 



- Li 2 (^^l-Li 2 (J-£±-)-U 2 {p.) 



P+ 



Li 2 



P- 



A+ -p. 
1 . o A + - p + 1 



+ Li 2 



A_ -p 4 



A_-p. 



logp_ log 



i - a_ ; ' 2 

l-p_ 1 



• i 2 ' = / / 4 ' i 2 P+ ~~ ^- ' 

+ - log — log 



1 



l-p + 



p+ 



- 2 P- 



P- 



2 /x - " 'I'- 1 1 P ' '\ 

log 



K - P- i , 

1 — p_ 2 1 — A_ 



, P- - 1 , A - 
log log ■ 



A_ — p~ A_ — p_ 



+ log 



A_ 



log 



-A. 



P- A_ - p_ 



7T 
+ "77" 



+ m 



2 log 



A+ - p -+log 1 -^ 



rpq—k 



K - P+ 

d d l 1 1 



1-P- 



(E.4) 



1 



(27r) d Z 2 (Z-p') 2 -m 2 (/+p) 2 -m 2 



iV(e) (m : 



1 I 1 



- <j - log J- + log — log - - log 2 £_ 



+ I log 2 e+ - log A' log L + Li 2 (-^j + Li 2 + I log 2 



it 

A' 



7T + Z7T 

6 

d d Z 1 



- + log 2 log A' 

e cr 3 



(E.5) 



2m 



1 



(27r) d Z 2 (p + /) 2 - m 2 (p + k + Z) 2 - m 2 
" ml0g ( 1 + S)} 



(E.6) 



d d l 1 



1 



(2ir) d / 2 (Z + /c) 2 (Z + p + /c) 2 -m 2 



= 7V(e) m — { — + -log— + log — log 1 + — 
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3m 



m 2 J 6 



- Li 2 ( -—77 ) - 77 ^ + ^ 



1 , / (7i \ m 2 
- + log 1 + ^ +21og — 



d d l 



r d, 
J CM 



1 



1 



1 



(2ir) d (l-p') 2 -m 2 (l+p-q) 2 -m 2 (I + p) 2 - m 2 



N(e) (m 2 )' 



W 7^-1 



- 1 



2«7T 



log 1-1 -log -1-1 



(E.7) 



(E.8) 



A partial check of the correctness of the above formulae can be performed in the 
following way. We consider first a check of B 2m - To this purpose, define I to be 



d d l 1 + A(l- k) 2 + B[{l+p-qf - m 2 } +C[(l + pf 



1 J^y 



m 



l 2 (l - k) 2 [(l+p- q) 2 - m 2 } [(/ + p) 2 - m 2 } 



I 



d d l 1 + B[q 2 -2p-q\+2l- [-Ak + B(p - q) + Cp] + l 2 (A + B + C) 



{2n) d l 2 (l - k) 2 [{l+p- q) 2 - m 2 ] [(I + p) 2 - m 2 ] 

If we impose that / has no infrared and collinear divergences, then 

' 1 + B[q 2 -2p-q] =0 

< 

^ k ■ [-Ak + B(p -q) + Cp] = 



Solving this system 



B = — 



C = 



0"2 
1 
0"! 



So I can be rewritten as 

1 



/ = B 2m + AT 2m 



-Ti 



0"2 



lm 



0"! 



lm 



where the primed quantities are the same as before, with the substitution p <-> p', 
that is a i <-> o" 2 . The integral / is now convergent and the cancellation of the di- 
vergent part of the right-hand side can be checked directly (both in the real part 
and in the absorptive one). As far as the finite terms are concerned, the integral / 
can be reduced to a one-variable integral, using Feynman parametrization, and then 
integrated numerically to check the identity. 
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With the same reasoning, we can check B 3m . We introduce the integral 

d d l 1 + A [{I + pf -m 2 ]+B[{l+p- qf -m 2 ] + C [(I - p'f - m 2 } 



r d l 



I 



{27i) d I 2 [{I + p) 2 - m 2 ] [{l + p- q) 2 - m 2 ] [{I - p') 2 - m 2 ] 

d d l 1 + B [q 2 -2p-q\+2l- [Ap + B{p - q) - Cp'\ + 1 2 {A + B + C) 



{2ix) d l 2 [{l + p) 2 -m 2 ][{l+p-q) 2 -m 2 ][{l-p') 2 -m 2 } 

This integral has only soft divergences, which can be removed if we require that 

1 



1 + B 



q 2 — 2p ■ q 



0^B = . 

0"2 



Thus / becomes 



I — B 3m + AT 2m T$ m k + CT 2 m . 
The rest of the check is the same as before. 



Appendix F: List of integrals for the soft contributions 

We now summarize the values of the integrals required to isolate the singular 
pieces of the four-jet cross section, in the soft-gluon limit 

.i_^_I logfc _ log ^ + i log2fc _^_ 2Li2 (_£)} + 0(£) 



I 2 (x, h) 



o dV Sl dvHl - v)rre ^+kv\ 2 



i-{-i-2io g (i + ^ + io g a +w 



±- r#(sin0)- 2e [ X dy f 1 dv[v{l-v)]-Y 
i\u Jo Jo Jo 



X 



y + hv y — c cos y/y \Jv + g : 



1 f 1 T T 



-Ix + h/x 



} + 0(e) 
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where I e , J^, I x and Iy x are finite quantities, defined by 



2 log sin j 

o (g — h) 2 + /i (ccos 4>) 



(2c(g + h)cos(f) c cos 

x 1 /7l 2 T T arctan //i 2 2 =^ + (9-h) log - 
[ v4(yf — cr cos z (p y/Ag — cr cos^ cp h t 



I x = I dt. l ° gt 



o (h + t) h t + g f _ c 2 t 



h/x = dt 



2 logo; — logt 7r 



' (l + ht) J(i + gt f _ c 2 t 



For the definition of the constants appearing in these integrals, see Section [4.4| and 
Appendix [E]. 



Appendix G: Results 

We implemented our analytical result in a FORTRAN program, which behaves 
like a "partonic" Monte Carlo generator, analogous to the program EVENT ||. We 
collect here some results obtained with our code, with which future users of the 
program may, eventually, compare their results. Furthermore, since for this kind of 
calculations it would be difficult to perform analytical comparisons, the only possible 
alternative is to choose a few shape variables, and compare numerical results, in the 
spirit of what has been done in ref. [JT7 for the case of the massless calculation. 



We include in these results only the contributions from cut graphs of A-type, that 
is to say, from cut graphs in which the weak current couples to the same heavy-flavour 
loop, and there is a single QQ pair in the final state, which is the really hard part 
of the calculation. For the contributions involving two heavy-quark pairs in the final 
state, it is easier to compare directly the value of the matrix elements squared (this 
part of our program was in fact checked in this way with the program of ref. ||10||). 
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We have chosen a set of shape variables for which it should be easy to obtain 
quite accurate numerical results. We have fixed the centre-of-mass energy to be 100 
GeV, and the mass of the heavy quark has been taken to be equal to 1, 10, 20 and 30 
GeV. We present separately the results for a hypothetical vector boson with purely 
axial or purely vector couplings, normalized to the massless total cross section at the 
zeroth order in a s . We have chosen the following shape variables: the thrust t, the 
c parameter, the mass of the heavy jet squared M% (according to the thrust axis), 
the energy-energy correlation EEC, the three-jet fractions according to the E, EM, 
JADE, and DURHAM schemes. For t, c, M| and EEC we present moments, instead 
of distributions, because they can be obtained with higher precision. For thrust, for 
example, we thus compute, according to the notation of Section [| 

TZ)^ V ~2tt ) Bty l A{ 



dTy/A (1 - t) n = ( A\r /A (n) + ) B v/A (n) . (G.l) 



We will further decompose 

By/A = B V / AtCA + B V i A Cf + B v / A Tf , (G.2) 

where the C A , Cf and Tp subscripts denote the CfC a , C f and fifCpTp colour com- 
ponents. For some shape variables, the presence of massive particles in the final state 
may introduce ambiguities in the definition, owing to the fact that, in the massless 
case, energy and momentum can be interchanged. We thus refer to the exact defini- 
tions given in ref. || for t, c, M| and in ref. |17| for the EEC. Moments are defined 
as 



dTv/Ae n = [^)A v/A {n)+{^) 2 B VIA {n) 



VIA ) = »v/aW f {^) ^ 

J dTy/A E ^f- cos fc % sin 2 +" % = (g) A v JZ(n, k) + (g) " B^(n, k) . 

where the sum runs over all the final particles. 

Clusters are defined in the following way. There is a resolution parameter y, which 
is computed for every pair of particles in the final state. One finds the pair for which 
y is minimum. If y < y cut the two particles are combined into a single pseudo-particle 
by adding up their four momenta. One thus computes 

dT v/A 5 Nx{ycutl3 = (J^J A$ /A (y cut ) + (^j B^ /A (y cut ) , (G.3) 
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where X stands for E, EM, JADE or DURHAM, and Nx{y C ut) is the number of pseudo- 
particles in the final state after the clustering procedure. The various clustering 
algorithms differ by the definition of the resolution parameter y 



E 



(Pi + Pj? 
q 2 



EM : 



JADE : 2 -VH1 -cos0«) 
q 2 



EiEq , 

-V i 

q 1 

DURHAM : 2 min (e?, E])-^(l - cos %) . (G.4) 

Observe that the E scheme is not infrared-safe if y cu t < m 2 /q 2 . In fact, in this case, 
the configuration made up of two heavy quarks plus a soft gluon cannot be reduced 
to two pseudo-particles, since the recombination parameter will fail the cut, for any 
pair containing a massive quark. The cancellation of soft divergences cannot therefore 
work for these values of the cut parameter. 

We have chosen the renormalization scale fi = E, and nf = 5. The results are 
given in Tables [I] to || 
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n 


m/£ = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 






1 


71.53 ±0.045 


57.14 ±0.024 




39.17± 0.013 


21.83 ±0.006 


2 


5.303 ± 0.005 


4.47 ±0.0034 




3.026 ±0.002 


1.524 ± 8 10~ 4 


3 


0.8056 ±0.0012 


0.6887 ± 8 10~ 4 




0.4582 ± 5 10~ 4 


0.2127 ± 2 10~ 4 


4 


0.1604 ± 3.4 10" 4 


0.1381 ± 2.4 10" 4 




0.09096 ± 1.4 10" 4 


0.04003 ± 6 10~ 5 


5 


0.03668 ± HO" 4 


0.03174 ± 7 10" 5 




0.02079 ± 4 10~ 5 


0.008845 ± 2 10" 5 






1 


-4.34 ±0.06 


-0.704 ±0.026 




2.449 ±0.01 


3.929 ±0.003 


2 


2.48 ±0.007 


1.78 ±0.004 




1.127 ± 0.0014 


0.699 ± 5 10~ 4 


3 


0.5348 ±0.002 


0.3883 ±0.001 




0.2316 ± 4 10~ 4 


0.1293 ± 1.2 10" 4 


4 


0.1241 ± 6 10~ 4 


0.09007 ± 2.7 10- 4 




0.0519 ± 1.1 10" 4 


0.02743 ± 3.4 10" 5 


5 


0.03116 ± 1.8 10" 4 


0.02255 ± 8 10" 5 




0.01261 ± 3.5 10" 5 


0.006392 ±110~ 5 




B'v.tM 


1 


-22.37 ±0.004 


-18.48 ±0.005 




-13.17 ±0.004 


-7.767 ±0.002 


2 


-1.552 ± 6 10~ 4 


-1.38±710~ 4 




-1.017±510~ 4 


-0.5844 ± 3 10~ 4 


3 


-0.2153 ± 1.3 10" 4 


-0.1958 ± 1.7 10" 4 




-0.1465 ± 1.2 10" 4 


-0.08272 ± 7 10" 5 


4 


-0.03864 ± 3.5 10- 5 


-0.03559 ± 4 10" 5 




-0.027 ± 3 10" 5 


-0.01515 ± 1.8 10" 5 


5 


-0.007897 ± 110- 5 


-0.007344 ± 1.3 10" 5 




-0.005647 ± 9 10- 6 


-0.003174 ± 5 10" 6 








cM 




1 


71.5 ±0.045 


54.27 ±0.023 




31.33 ±0.01 


12.46 ±0.0035 


2 


5.301 ±0.005 


4.284 ±0.003 




2.501 ±0.0015 


0.9336 ± 5 10~ 4 


3 


0.8054 ±0.0012 


0.6637 ± 8 10~ 4 




0.3874 ± 4 10~ 4 


0.138 ± 1.2 10~ 4 


4 


0.1603 ± 3.4 10" 4 


0.1336 ± 2.3 10" 4 




0.07807 ± 110~ 4 


0.02708 ± 3.6 10" 5 


5 


0.03667 ± HO" 4 


0.03077 ± 7 10" 5 




0.01802 ± 3.5 10" 5 


0.006162 ± 1.1 10" 5 






B X 


cA n ) 




1 


-4.3 ±0.06 


0.5 ±0.025 




3.942 ±0.008 


3.722 ±0.002 


2 


2.481 ±0.007 


1.801 ±0.0036 




1.109 ±0.0012 


0.5694 ± 3 10~ 4 


3 


0.535 ±0.002 


0.3871 ± 9 10~ 4 




0.221 ± 3 10~ 4 


0.1041 ± 8 10~ 5 


4 


0.1241 ± 6 10~ 4 


0.0893 ± 2.6 10" 4 




0.04897 ± 9 10~ 5 


0.02206 ± 2 10~ 5 


5 


0.03116 ± 1.8 10" 4 


0.02228 ± 8 10" 5 




0.01181 ± 3 10~ 5 


0.005136 ± 7 10" 6 








T F ( n ) 




1 


-22.36 ±0.004 


-17.5 ±0.005 




-10.49 ±0.003 


-4.457 ±0.0013 


2 


-1.552 ± 6 10~ 4 


-1.311±710~ 4 




-0.8247 ± 4 10~ 4 


-0.3533 ± 1.7 10- 4 


3 


-0.2152 ± 1.3 10" 4 


-0.1861 ± 1.6 10" 4 




-0.1197±910~ 5 


-0.05156 ± 4 10" 5 


4 


-0.03862 ± 3.5 10- 5 


-0.03382 ± 4 10" 5 




-0.0221 ± 2.5 10" 5 


-0.009594 ± 1 10~ 5 


5 


-0.007893 ± 110- 5 


-0.00697 ± 1.2 10" 5 




-0.004616 ± 7 10" 6 


-0.002024 ± 3 10" 6 



Table 1: The thrust t. 
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n 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 






D V,C A 


(») 




1 


68.32 ±0.05 


56.71 ±0.03 




38.45 ±0.02 


19.46 ±0.008 


2 


4.591 ±0.006 


4.014 ±0.004 




2.642 ± 0.0035 


1.029 ±0.0024 


3 


0.6336 ±0.0013 


0.5549 ±0.001 




0.3565 ±0.0013 


0.1072 ±0.001 


4 


0.1156 ± 3.6 10 -4 


0.1002 ± 3 10~ 4 




0.0641 ± 5 10~ 4 


0.0167 ± 4.5 10~ 4 


5 


0.02447 ± 1 10~ 4 


0.02084 ±110~ 4 




0.01346 ± 1.8 10" 4 


0.0037 ± 2 10~ 4 






R M h 


(n) 




1 


-21.39 ±0.05 


-10.55 ±0.026 




-1.143 ±0.01 


3.166 ±0.004 


2 


0.218 ±0.007 


0.4605 ±0.004 




0.7282 ±0.002 


0.736 ±0.0011 


3 


0.0784 ±0.002 


0.119± 0.001 




0.162 ± 6 10~ 4 


0.1834 ± 5 10~ 4 


4 


0.0165 ± 6 10~ 4 


0.02626 ± 3 10~ 4 




0.03773 ± 2 10~ 4 


0.055 ± 2 10~ 4 


5 


0.00365 ± 1.8 10~ 4 


0.00619 ±110~ 4 




0.00959 ± 8 10~ 5 


0.01829 ± 9 10~ 5 






B Mh 

V,T F 


(n) 




1 


-23.2 ±0.005 


-20.04 ±0.006 




-14.32 ±0.004 


-7.809 ±0.0025 


2 


-1.704 ± 6 10~ 4 


-1.563 ± 8 10~ 4 




-1.153±710~ 4 


-0.5961 ± 5 10~ 4 


3 


-0.2505 ± 1.4 10" 4 


-0.2308 ± 2 10~ 4 




-0.1753 ± 2 10~ 4 


-0.1014 ± 2 10~ 4 


4 


-0.04758 ± 3.6 10 -5 


-0.04375 ± 5 10~ 5 




-0.03474 ± 6 10" 5 


-0.02601 ± 8 10~ 5 


5 


-0.0103 ± 1 10~ 5 


-0.009445 ± 1.5 10" 5 




-0.007984 ± 2 10" 5 


-0.0082 ± 3 10~ 5 






B A,C A 


(n) 




1 


68.29 ±0.05 


53.84 ±0.03 




30.76 ±0.016 


10.99 ±0.006 


2 


4.59 ±0.006 


3.841 ±0.004 




2.187 ±0.004 


0.6034 ±0.0018 


3 


0.6334 ±0.0013 


0.533 ±0.001 




0.3026 ±0.0015 


0.0638 ± 8 10~ 4 


4 


0.1156 ± 3.6 10~ 4 


0.0965 ± 3 10~ 4 




0.0553 ± 6 10~ 4 


0.0102 ± 3.4 10~ 4 


5 


0.02447 ± 1 10~ 4 


0.02013 ±110~ 4 




0.01173 ± 2 10~ 4 


0.00246 ± 1.4 10~ 4 






r,M h 

A, Cf 


(n) 




1 


-21.35 ±0.05 


— 8.804 ± 0.025 




1.276 ±0.009 


3.617 ±0.0026 


2 


0.22 ±0.007 


0.546 ±0.004 




0.8205 ±0.0017 


0.7156 ± 7 10~ 4 


3 


0.0787 ±0.002 


0.1304 ±0.001 




0.1712 ± 6 10~ 4 


0.1817±310~ 4 


4 


0.0166 ± 6 10~ 4 


0.02845 ± 3 10~ 4 




0.03897 ± 2 10~ 4 


0.05504 ± 1.2 10" 4 


5 


0.00366 ± 1.8 10- 4 


0.0067 ± 9 10~ 5 




0.00974 ± 8 10~ 5 


0.01825 ± 5 10" 5 






A, Tf 


(n) 




1 


-23.19 ±0.005 


— 19.01 ± 0.005 




-11.52 ±0.0033 


-4.617 ±0.0015 


2 


-1.704 ± 6 10~ 4 


— 1.491 ± 8 10 4 




-0.9601 ± 6 10~ 4 


-0.4038 ± 3.6 10- 4 


3 


-0.2504 ± 1.4 10" 4 


-0.2207 ± 1.7 10" 4 




-0.149 ± 1.5 10" 4 


-0.07702 ± 1.4 10" 4 


4 


-0.04756 ± 3.6 10" 5 


-0.04188 ± 4.5 10" 5 




-0.02987 ± 5 10- 5 


-0.02099 ± 5 10" 5 


5 


-0.0103 ± 1 10" 5 


-0.009048 ± 1.3 10~ 5 




-0.006893 ± 1.6 10" 5 


-0.00675 ± 2 10" 5 



Table 2: The mass of the heavy jet squared M%. 
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n 


m/£ = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 




B c v ,c A (n) 


1 


303.3 ±0.18 


233.9 ±0.09 


155.2 ±0.05 


84.7 ±0.023 


2 


70.97 ±0.06 


58.25 ±0.03 


38.33 ±0.018 


19.18 ±0.007 


3 


30.55 ±0.03 


25.56 ±0.02 


16.63 ±0.01 


7.789 ±0.004 


4 


16.36 ±0.02 


13.85 ±0.012 


8.955 ±0.007 


4.008 ± 0.003 


5 


9.784 ±0.012 


8.341 ±0.008 


5.379 ±0.005 


2.335 ±0.002 




B c v ,cM 


1 


-36.2 ±0.2 


-13.58 ±0.1 


4.84 ± 0.04 


13.41 ±0.013 


2 


30.42 ±0.07 


20.29 ±0.036 


12.49 ±0.013 


7.845 ±0.004 


3 


20.1 ±0.04 


13.69 ±0.02 


7.741 ± 0.008 


4.23 ±0.003 


4 


13.19 ±0.025 


9.04 ±0.013 


4.892 ± 0.005 


2.479 ±0.002 


5 


9.024 ±0.017 


6.185 ±0.009 


3.232 ±0.004 


1.543 ±0.0013 




B c v .tM 


1 


-95.12 ±0.016 


-75.92 ±0.02 


-52.26 ±0.015 


-30.07 ±0.009 


2 


-20.81 ±0.006 


-18.02 ±0.007 


-12.82 ±0.005 


-7.206 ±0.003 


3 


-8. 117 ±0.003 


-7.227 ±0.004 


-5.243 ±0.003 


-2.915 ±0.0018 


4 


-3.859 ±0.002 


-3.495 ±0.0027 


-2.581 ±0.002 


-1.435 ±0.0012 


5 


-2. 007 ±0.0013 


-1.842 ±0.0018 


-1.387 ±0.0013 


-0.7771 ± 8 10~ 4 






B C A 






1 


303.2 ±0.18 


222 ± 0.09 


123.9 ±0.04 


48.14 ±0.014 


2 


70.95 ±0.06 


55.73 ±0.03 


31.48 ±0.014 


11.58 ±0.0045 


3 


30.54 ±0.03 


24.59 ±0.018 


13.95 ±0.008 


4.942 ± 0.0027 


4 


16.36 ±0.02 


13.36 ±0.01 


7.62 ±0.005 


2.643 ±0.0018 


5 


9.782 ±0.012 


8.071 ±0.007 


4.628 ±0.004 


1.587 ±0.0013 






B A 


,Cf(") 




1 


-36.04 ±0.2 


-8.24 ±0.1 


11.47 ±0.03 


13.16 ±0.008 


2 


30.44 ±0.07 


20.63 ±0.034 


12.45 ±0.011 


6.425 ± 0.0027 


3 


20.11 ±0.04 


13.64 ±0.02 


7.389 ±0.007 


3.394 ±0.0017 


4 


13.19 ±0.025 


8.94 ±0.013 


4.589 ±0.004 


1.975 ±0.001 


5 


9.024 ±0.017 


6.091 ±0.009 


2.998 ±0.003 


1.223 ± 8 10~ 4 






a A 






1 


-95.07 ±0.016 


-71.88 ±0.02 


-41.59 ±0.012 


-17.2 ±0.005 


2 


-20.8 ±0.006 


-17.11 ±0.007 


-10.36 ±0.004 


-4.312 ±0.0017 


3 


-8.113 ±0.003 


-6.865 ±0.004 


-4.263 ± 0.0025 


-1.788 ±0.001 


4 


-3.857 ±0.002 


-3.317 ±0.0025 


-2.099 ±0.0016 


-0.8903 ± 7 10~ 4 


5 


-2.006 ±0.0013 


-1.745 ±0.0017 


-1.124 ±0.001 


-0.4836 ± 4.5 10" 4 



Table 3: The c parameter. 
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n 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 




£y E c>>0) 


U 


202.3 ±0.12 


154.4 ±0.06 


97.08 ±0.03 


46.81 ±0.013 


i 
1 


143.1 ±0.1 


117 ±0.05 


77.18 ±0.026 


38.2 ±0.01 


Z 


115.9 ±0.09 


97.55 ±0.04 


66.33 ±0.022 


33.55 ±0.009 


Q 
O 


99.79 ±0.08 


85.27 ±0.04 


59.11 ±0.02 


30.38 ±0.008 


A 
ft 


88.92 ±0.08 


76.65 ±0.04 


53.84 ±0.02 


28.01 ±0.007 


5 


80.95 ±0.07 


70.18 ±0.04 


49.76 ±0.018 


26.13 ±0.007 




^ E c c F (n,0) 


u 


-24.13 ±0.14 


-8.83 ±0.07 


3.567 ±0.023 


8.224 ±0.007 


i 
i 


-9.82 ±0.12 


-3.46 ±0.05 


3.814 ±0.02 


6.861 ±0.006 


o 
z 


-5.96 ±0.1 


-1.79 ±0.05 


3.55 ±0.016 


5.99 ±0.005 


Q 
O 


-4.33 ±0.1 


-1.06 ±0.04 


3.284 ±0.015 


5.387 ±0.005 


/I 


-3.46 ±0.1 


-0.664 ±0.04 


3.06 ±0.014 


4.937 ±0.0045 


5 


-2.92 ±0.1 


-0.43 ±0.04 


2.873 ±0.013 


4.585 ± 0.004 






u 


-63.43 ±0.01 


-50.32 ± 0.014 


-33.15 ±0.009 


-17.12 ±0.005 


i 
l 


-45.25 ±0.009 


-38.32 ±0.01 


-26.45 ±0.007 


-13.98 ±0.004 


o 
Z 


-36.83 ±0.008 


-32.05 ±0.01 


-22.76 ±0.006 


-12.27 ±0.0035 


q 
o 


-31.8 ±0.008 


-28.07 ±0.009 


-20.3 ±0.006 


-11.11 ±0.003 


/l 
ft 


-28.39 ±0.007 


-25.26 ±0.008 


-18.5 ±0.005 


-10.23 ±0.003 


5 


-25.88 ±0.007 


-23.15 ±0.008 


-17.1 ±0.005 


-9.544 ± 0.0027 






U 


202.2 ±0.12 


146.6 ±0.06 


77.59 ±0.027 


26.66 ±0.008 


1 


143 ±0.1 


111 ±0.05 


61.68 ±0.02 


21.72 ±0.006 


z 


115.8 ±0.09 


92.62 ±0.04 


52.98 ±0.018 


19.05 ±0.006 


q 
o 


99.75 ±0.08 


80.96 ±0.04 


47.2 ±0.016 


17.23 ±0.005 


/i 
ft 


88.87 ±0.08 


72.77 ±0.04 


42.97 ±0.015 


15.87 ±0.005 


5 


80.91 ±0.07 


66.62 ±0.035 


39.71 ±0.014 


14.79 ±0.005 




^ E c C >,0) 


n 
u 


-24.04 ±0.14 


-5.3 ±0.06 


7.765 ± 0.02 


8.116 ±0.0045 


i 
i 


-9.76 ±0.12 


-0.89 ±0.05 


7.004 ±0.016 


6.709 ±0.004 


o 


-5.9 ±0.1 


0.32 ±0.045 


6.252 ±0.014 


5.85 ±0.003 


Q 

o 


-4.28 ±0.1 


0.78 ±0.04 


5.676 ±0.013 


5.259 ±0.003 


4 


-3.41 ±0.1 


0.98 ±0.04 


5.23 ±0.012 


4.82 ±0.003 


5 


-2.88 ±0.1 


1.07 ±0.04 


4.872 ±0.011 


4.476 ±0.003 




B™g(n,0) 





-63.4 ±0.01 


-47.66 ±0.013 


-26.48 ±0.008 


-9.928 ±0.003 


1 


-45.22 ±0.009 


-36.29 ±0.01 


-21.1 ±0.006 


-8.091 ± 0.0024 


2 


-36.81 ±0.008 


-30.35 ±0.009 


-18.15 ±0.005 


-7.085 ±0.002 


3 


-31.79 ±0.008 


-26.57 ±0.008 


-16.18 ±0.005 


-6.401 ±0.002 


4 


-28.38 ±0.007 


-23.91 ±0.008 


-14.73 ±0.004 


-5.889 ±0.0017 


5 


-25.87 ±0.007 


-21.91 ±0.007 


-13.61 ±0.004 


-5.486 ±0.0016 



Tahlp A- Tlnp pnprcfv- pnprcrv rnrrplatinn h = D 
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n 


m/£ = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 




^ A (n,l) 


A 
U 


-26.06 ±0.025 


-21.45 ± 0.014 


-14.12 ±0.008 


-7.006 ±0.003 


1 
1 


-9.56 ±0.016 


-8. 967 ±0.01 


-6.413 ±0.006 


-3.152 ± 0.0023 


Z 


-4.827± 0.013 


-4.892 ±0.008 


-3.785 ±0.0045 


-1.898 ±0.002 


q 
O 


-2.894 ±0.011 


-3.086 ±0.008 


-2.549 ±0.004 


-1.316 ±0.0017 


A 


-1.93 ±0.01 


-2.133 ±0.007 


-1.858 ±0.0036 


-0.9876 ±0.0015 


5 


-1.38 ±0.01 


-1.57 ±0.007 


-1.427 ±0.0034 


-0.779 ±0.0014 




B$%(n,l) 


A 
U 


8.73 ±0.03 


3.607 ±0.015 


-0.625 ±0.006 


-2.075 ±0.002 


1 
1 


0.593 ±0.02 


0.272 ±0.01 


-0.944 ±0.004 


-1.345 ±0.0014 


Z 


-0.588 ±0.02 


-0.327 ±0.009 


-0.786 ±0.0034 


-0.9457 ±0.0012 


Q 

o 


-0.76 ±0.018 


-0.441 ±0.008 


-0.6335 ±0.003 


-0.7102 ±0.001 


4: 


-0.73 ±0.017 


-0.437 ±0.007 


-0.5177 ±0.003 


-0.5588 ± 9 10~ 4 


5 


-0.66 ±0.016 


-0.403 ±0.007 


-0.431 ±0.0026 


-0.4546 ± 8 10~ 4 






A 
U 


10.75 ±0.0026 


9.205 ±0.003 


6.342 ±0.0023 


3.381 ±0.0013 


1 
1 


4.47 ± 0.0018 


4.256 ± 0.0023 


3.179 ±0.0016 


1.718 ± 9 10~ 4 


Z 


2.5 ±0.0014 


2.506 ±0.002 


1.994 ±0.0013 


1.108 ± 7 10~ 4 


Q 

o 


1.629 ±0.0012 


1.679 ±0.0016 


1.4 ±0.001 


0.7997 ± 6 10~ 4 


/I 


1.162 ±0.001 


1.218 ±0.0014 


1.052 ± 9 10~ 4 


0.616 ± 5 10~ 4 


5 


0.8805 ±0.001 


0.933 ±0.0013 


0.8273 ± 8 10~ 4 


0.4948 ± 5 10~ 4 






A 
U 


-26.05 ±0.025 


-20.4 ±0.014 


-11.39 ±0.006 


-4.12 ±0.002 


1 
1 


-9.555 ±0.016 


-8.558 ±0.01 


-5.234 ±0.004 


-1.907 ± 0.0014 


z 


-4.825 ±0.013 


-4.68 ±0.008 


-3.108 ±0.0036 


-1.165± 0.0012 


Q 
O 


-2.894 ±0.011 


-2.957 ±0.007 


-2.102 ±0.003 


-0.8135 ±0.001 




-1.93 ±0.01 


-2.047 ±0.007 


-1.535 ±0.003 


-0.6129 ± 9 10~ 4 


5 


-1.38 ±0.01 


-1.508 ±0.006 


-1.182 ± 0.003 


-0.4846 ± 9 10~ 4 




^ E c C >,l) 


A 
U 


8.71 ±0.03 


2.855 ±0.014 


-1.525 ±0.005 


-1.98 ±0.0012 


1 
1 


0.586 ±0.02 


-0.025 ±0.01 


-1.315 ±0.0034 


-1.228 ± 9 10~ 4 


o 

z 


-0.59 ±0.02 


-0.486 ±0.008 


-0.997 ±0.003 


-0.853 ± 7 10~ 4 


Q 

o 


-0.763 ±0.018 


-0.541 ±0.007 


-0.7733 ±0.0025 


-0.6373 ± 6 10~ 4 


4 


-0.732 ±0.017 


-0.506 ±0.007 


-0.6186 ±0.0023 


-0.4999 ± 6 10~ 4 


5 


-0.66 ±0.016 


-0.454 ±0.007 


-0.508 ±0.002 


-0.4059 ± 5 10~ 4 









10.75 ± 0.0026 


8.765 ±0.003 


5.171 ±0.002 


2.054 ± 7 10~ 4 


1 


4.468 ±0.0018 


4.06 ±0.002 


2.614 ±0.0013 


1.07 ± 5 10~ 4 


2 


2.499 ±0.0014 


2.393 ±0.0018 


1.646 ±0.001 


0.6961 ± 4 10~ 4 


3 


1.628 ±0.0012 


1.603 ±0.0015 


1.157±810~ 4 


0.5043 ± 3.6 10- 4 


4 


1.162 ±0.001 


1.164 ±0.0013 


0.8701 ± 7 10~ 4 


0.3888 ± 3 10~ 4 


5 


0.88 ±0.001 


0.8916 ±0.0012 


0.6845 ± 7 10~ 4 


0.3123 ± 3 10~ 4 



Tahlp Jv Tlip pnprcrv— pnprcrv rnrrplatinn FYR(~! h = 1 
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Vcut 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 




B v,c A (V cut ) 


0.01 


1357 ± 1.5 








0.05 


341.3 ±0.5 


374 ±0.3 


817.1 ±0.3 




0.10 


133.7 ±0.3 


138.6 ±0.17 


164.3 ±0.1 


459.5 ±0.17 


0.15 


60.75 ±0.18 


63.1 ±0.13 


64.29 ±0.07 


79.06 ±0.04 


0.20 


27.13 ±0.16 


29 ±0.12 


29 ±0.05 


24.54 ±0.027 




B V,C F iVcut) 


0.01 


-88 ±4 








0.05 


129.8 ± 1.3 


249.1 ±0.36 


0.84 ±0.4 




0.10 


71.8 ±0.7 


120.1 ±0.2 


94.69 ±0.07 


47.95 ±0.13 


0.15 


37.5 ±0.5 


53.8 ±0.13 


59.02 ±0.045 


41.25 ±0.025 


0.20 


18.1 ±0.34 


24.26 ±0.13 


32.27 ±0.03 


21.58 ±0.014 




B V,T F (.Vcut) 


0.01 


-452.3 ±0.13 








0.05 


-103.8 ±0.04 


-121.3 ±0.05 


-268.6 ±0.12 




0.10 


-38.34 ±0.02 


-42.14 ±0.03 


-54.59 ±0.022 


-154.4 ±0.07 


0.15 


-16.56 ±0.015 


-17.79 ±0.02 


-21.16 ±0.015 


-28. 17 ± 0.012 


0.20 


-6.871 ±0.009 


-7.424 ±0.013 


-8.984 ±0.012 


-9.268 ±0.008 






0.01 


1356 ± 1.5 








0.05 


341.2 ±0.5 


355.5 ±0.3 


643 ±0.27 




0.10 


133.6 ±0.3 


132.8 ±0.17 


133.1 ±0.09 


254.5 ±0.09 


0.15 


60.9 ±0.24 


60.83 ±0.13 


53.63 ±0.06 


46.55 ±0.03 


0.20 


27.1 ±0.14 


27.83 ±0.08 


25.03 ±0.05 


15.4 ±0.016 




b a,c f (y<™t) 


0.01 


-88 ±3.6 








0.05 


130 ± 1.3 


243.1 ±0.3 


32.1 ±0.3 




0.10 


71.8 ±0.7 


117.2 ±0.2 


85.79 ±0.06 


49.67 ±0.07 


0.15 


37.45 ±0.5 


52.54±0.12 


52.99 ±0.04 


31.7 ±0.016 


0.20 


18.1 ±0.33 


23.8 ±0.09 


29.08 ±0.03 


17.24 ±0.01 




B A,T F (ycut) 


0.01 


-452.1 ±0.13 








0.05 


-103.8 ±0.04 


-115.4 ±0.045 


-212.6 ±0.1 




0.10 


-38.32 ±0.02 


-40.23 ±0.027 


-44.67 ±0.02 


-86.57 ±0.04 


0.15 


-16.55 ±0.015 


-17.01 ±0.02 


-17.8 ±0.013 


-17.28 ±0.008 


0.20 


-6.868 ±0.009 


-7.098 ±0.012 


-7.747 ±0.01 


-6.237 ±0.005 



Table 6: The E clustering algorithm. 
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Vent 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 






0.01 


1340 ± 1.7 


1135 ±0.8 


803.4 ±0.4 


455. 7± 0.18 


0.05 


333.2 ±0.5 


291.2 ±0.4 


203.4 ±0.18 


101.6 ±0.07 


0.10 


126.3 ±0.3 


111.8 ±0.23 


73.76 ±0.11 


29.24 ±0.04 


0.15 


54.24 ±0.3 


48 ±0.17 


30.47 ±0.12 


8.176 ±0.027 


0.20 


22.02 ±0.2 


18.85 ±0.14 


10.77 ±0.08 


1.52 ±0.024 




B v!c F iVcut) 


0.01 


-327 ±4 


-193.4 ± 1.5 


-58.2 ±0.4 


37.07 ±0.14 


0.05 


74.5 ± 1.4 


75.8 ±0.6 


68.44 ±0.14 


43.09 ±0.04 


0.10 


39.6 ±0.8 


40.85 ±0.3 


35.2 ±0.16 


20.79 ±0.04 


0.15 


18.5 ±0.5 


18.47 ±0.27 


15.9 ±0.06 


9.353 ±0.015 


0.20 


7.46 ±0.3 


6.45 ±0.5 


5.98 ±0.05 


3.57 ±0.011 






0.01 


-453.4 ±0.14 


-391.2 ±0.13 


-276.1 ±0.12 


-157±0.07 


0.05 


-106.9 ±0.04 


-98.04 ±0.05 


-70.61 ±0.03 


-36.9 ±0.015 


0.10 


-40.71 ±0.025 


-37.79 ±0.03 


-27.29 ±0.018 


-12.24 ±0.01 


0.15 


-18.25 ±0.02 


-16. 7± 0.022 


-11.99 ±0.015 


-4.578 ±0.006 


0.20 


-7.97 ±0.01 


-7.125 ±0.013 


-4.687 ±0.01 


-1.53 ±0.005 




B A M c A (ycut) 


0.01 


1339 ± 1.7 


1072 ±0.9 


629.5 ±0.3 


250.8 ±0.14 


0.05 


332.8 ±0.5 


277.2 ±0.35 


163.4 ±0.17 


58.88 ±0.07 


0.10 


126.3 ±0.3 


107.1 ±0.26 


61.04 ±0.16 


17.97 ±0.03 


0.15 


54.2 ±0.2 


46.08 ±0.2 


25.54 ±0.1 


5.185 ±0.02 


0.20 


21.97 ±0.2 


17.9 ±0.25 


9.23 ±0.08 


1.013 ±0.014 




B A% F (Vcut) 


0.01 


-326.5 ±4 


-164.4 ± 1.2 


-15.6 ±0.3 


43.13 ±0.07 


0.05 


74.6 ± 1.4 


78.6 ±0.5 


66.16 ±0.13 


33.78 ±0.022 


0.10 


39.5 ±0.8 


40.8 ±0.6 


34.34 ±0.08 


17.18 ±0.013 


0.15 


18.4 ±0.5 


18.7 ±0.2 


15.7 ±0.07 


8.532 ±0.012 


0.20 


7.36 ±0.35 


7.29 ±0.14 


5.96 ±0.05 


3.529 ±0.007 




B A^r F (y^) 


0.01 


-453.1 ±0.14 


-370. 7± 0.13 


-219.1 ±0.1 


-88.5 ±0.04 


0.05 


-106.8 ±0.04 


-93.35 ±0.04 


-57.47 ±0.023 


-22.27 ±0.009 


0.10 


-40.69 ±0.024 


-36.04 ±0.03 


-22.83 ±0.016 


-8.066 ±0.006 


0.15 


-18.23 ±0.017 


-15.98 ±0.02 


-10.21 ±0.01 


-3.347 ±0.004 


0.20 


-7.965 ±0.01 


-6.83 ±0.013 


-4.033 ±0.008 


-1.247 ±0.0035 



Table 7: The EM clustering algorithm. 
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Vent 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 




B v A Sf(ycut) 


0.01 


1352 ±2 


1088 ± 1.6 


727.2 ±0.5 


389.2 ±0.26 


0.05 


328.6 ±0.6 


275 ±0.5 


180 ±0.2 


82.34 ±0.09 


0.10 


124.2 ±0.3 


106 ±0.26 


66.06 ±0.15 


24.91 ±0.05 


0.15 


54.1 ±0.3 


46.84 ±0.26 


28.77 ±0.1 


8.675 ±0.03 


0.20 


22.86 ±0.2 


19.84 ±0.14 


12.34 ±0.09 


2.98 ±0.03 




R JADE/„, \ 


0.01 


-787 ±5 


-423.2 ± 1.5 


-112.4 ±0.4 


24.27 ± 0.14 


0.05 


11.7 ± 1.4 


25.2 ±0.7 


35.2 ±0.4 


29.65 ±0.04 


0.10 


27.8 ±0.8 


29.25 ±0.33 


25.1 ±0.09 


15.33 ±0.027 


0.15 


17.45 ±0.5 


16.64 ±0.24 


13.61 ±0.08 


7.885 ±0.018 


0.20 


8.7 ±0.3 


7.7 ±0.14 


6.6 ±0.06 


3.796 ±0.012 






0.01 


-470 ±0.16 


-380.9 ±0.16 


-255.8 ±0.11 


-138.5 ±0.06 


0.05 


-110.8 ±0.05 


-96.21 ±0.05 


-65.03 ±0.03 


-32.06 ±0.015 


0.10 


-41.68 ±0.025 


-37.65 ±0.03 


-25.47 ±0.02 


-11.21 ±0.01 


0.15 


-18.34 ±0.017 


-16.84 ±0.025 


-11.73 ±0.013 


-4.688 ±0.007 


0.20 


-7.834 ±0.01 


-7.252 ±0.013 


-5.176 ±0.01 


-1.955 ±0.005 






0.01 


1350 ±2 


1027 ± 1 


569.9 ±0.5 


214.1 ±0.3 


0.05 


328.6 ±0.6 


261.9 ±0.4 


144.4 ±0.17 


47.69 ±0.05 


0.10 


124.1 ±0.3 


101.6 ±0.3 


54.78 ±0.1 


15.35 ±0.03 


0.15 


54.15 ±0.27 


45.25 ±0.2 


24.22 ±0.09 


5.72 ±0.02 


0.20 


23 ±0.2 


19.12 ±0.16 


10.4 ±0.12 


2.113 ±0.02 




Bi A S*(ycut) 


0.01 


-785 ± 5 


-384 ± 1.8 


-59.6 ±0.3 


33.72 ±0.1 


0.05 


11.7 ± 1.4 


30.6 ±0.5 


38.6 ±0.2 


24.55 ±0.04 


0.10 


27.9 ±0.8 


30 ±0.36 


25.05 ±0.11 


12.84 ±0.02 


0.15 


17.3 ±0.5 


17.35 ±0.23 


13.73 ±0.06 


6.91 ±0.01 


0.20 


8.7 ±0.3 


8.05 ±0.17 


6.66 ±0.045 


3.461 ±0.01 




Bi A ^(y C ut) 


0.01 


-469.7 ±0.16 


-360.5 ±0.2 


-202.8 ±0.09 


-78.07 ±0.036 


0.05 


-110.8 ±0.05 


-91.5 ±0.045 


-52.62 ±0.03 


-19.17 ±0.009 


0.10 


-41.66 ±0.025 


-35.93 ±0.03 


-21.09 ±0.015 


-7.253 ±0.006 


0.15 


-18.33 ±0.017 


-16.1 ±0.02 


-9.935 ±0.01 


-3.29 ±0.004 


0.20 


-7.826 ±0.01 


-6.96 ±0.014 


-4.448 ±0.008 


-1.476 ±0.0035 



Table 8: The JADE clustering algorithm. 
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Vent 


m/E = 0.01 


m/E = 0.1 


m/E = 0.2 


m/E = 0.3 






0.01 


443 ± 1.6 


329.9 ±0.5 


181.6 ±0.2 


71.09 ±0.07 


0.05 


116 ±0.3 


97.26 ±0.22 


52.4 ±0.12 


12.98 ±0.05 


0.10 


45.25 ±0.2 


39.53 ±0.2 


22.5 ±0.09 


3.62 ±0.03 


0.15 


20.4 ±0.18 


17.82 ±0.13 


10.37 ±0.08 


1.16 ±0.026 


0.20 


8.7 ±0.12 


7.61 ±0.1 


4.63 ±0.05 


0.41 ±0.02 




Bv^CpiVcut) 


0.01 


-128.3 ±2 


-51.2 ± 0.6 


10.43 ±0.18 


23.01 ±0.04 


0.05 


3.8 ± 1 


11.9 ±0.4 


16.8 ±0.1 


10.58 ±0.024 


0.10 


6.1 ±0.5 


8.84 ±0.2 


9.98 ±0.08 


5.444 ±0.017 


0.15 


4 ±0.3 


4.6 ±0.15 


5.32 ±0.05 


2.824 ±0.01 


0.20 


1.44 ±0.24 


2.09 ±0.1 


2.44 ±0.03 


1.039 ±0.009 




Buyout) 


0.01 


-162.5 ±0.07 


-127.3 ±0.06 


-73.77 ±0.03 


-31.63 ±0.014 


0.05 


-42.66 ±0.026 


-37.61 ±0.034 


-23.03 ±0.02 


-7.768 ±0.008 


0.10 


-17.37 ±0.015 


-15.8 ±0.02 


-10.49 ±0.014 


-2.996 ±0.006 


0.15 


-8.024 ±0.011 


-7.353 ±0.014 


-5.046 ±0.01 


-1.302 ±0.005 


0.20 


-3.503 ±0.008 


-3.23 ±0.01 


-2.206 ±0.008 


-0.427 ±0.004 




BTcJvcut) 


0.01 


443.6 ± 1 


313 ±0.6 


146.1 ±0.17 


41.43 ±0.06 


0.05 


115.4 ±0.6 


93 ±0.3 


43.97 ±0.12 


8.37 ±0.024 


0.10 


45.07 ±0.22 


37.75 ±0.18 


19.15 ± 0.1 


2.554 ±0.02 


0.15 


20.48 ±0.16 


17.2 ±0.14 


9.11 ±0.06 


0.887 ±0.016 


0.20 


8.77±0.12 


7.3 ±0.1 


4.04 ±0.05 


0.34 ±0.018 




BT$ F {ycut) 


0.01 


-127 ±2 


-42.3 ±0.7 


20 ±0.16 


21.21 ±0.03 


0.05 


3.9 ±0.9 


14.5 ±0.35 


18.2 ±0.1 


9.74 ±0.017 


0.10 


6 ±0.5 


9.77 ±0.2 


10.44 ±0.07 


5.213 ±0.011 


0.15 


4.04 ±0.3 


4.84 ±0.2 


5.51 ±0.06 


2.768 ±0.007 


0.20 


1.47 ±0.25 


2.23 ±0.1 


2.5 ±0.036 


1.012 ±0.006 




B^Mut) 


0.01 


-162.4 ±0.06 


-121.3 ±0.05 


-60.29 ±0.025 


-19.36 ±0.009 


0.05 


-42.65 ±0.026 


-36 ±0.03 


-19.5 ±0.016 


-5.433 ±0.005 


0.10 


-17.36 ±0.015 


-15.15 ±0.02 


-9.073 ±0.014 


-2.345 ±0.004 


0.15 


-8.026 ±0.012 


-7.06 ±0.013 


-4.407 ±0.009 


-1.097 ±0.0034 


0.20 


-3.499 ±0.008 


-3.093 ±0.008 


-1.934 ±0.006 


-0.369 ± 0.0025 



Table 9: The DURHAM clustering algorithm. 
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